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SUMMARY 
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The NASA Aircraft Noise Prediction Program (ANOPP) Propeller Analysis 
System is a set of computational modules for predicting the aerodynamics, 
performance, and noise of propellers. Propeller blade geometry is given in 
terms of blade surface coordinates derived from a Joukowski transform of the 
blade sections. Potential flow around the blade sections is computed by 
Theodorsen's method using the Kutta condition to fix the circulation. Blade 
boundary layers are computed using the Holstein-Bohlen method in the laminar 
region and the Truckenbrodt method in the turbulent region. Profile drag is 
predicted by the method of Young and Squires. Performance and induced flow 
are computed by Lock's method with the Prandtl circulation function near the 
blade tip. Discrete tone noise is predicted from blade shape and aerodynamic 
loads using Farassat's methods: the blade surface integral method for 

subsonic propellers and the collapsing sphere method for transonic 
propellers. Broadband trailing edge noise is computed by Schl inker and 
Amiet's method. The results of this prediction system are compared to 
measurements on two propellers: one subsonic and one transonic. Nearfield 

levels on the subsonic propeller are accurately predicted if the predicted 
power coefficient is adjusted to match the measured power coefficient of the 
propeller. The lower frequency harmonics of the subsonic propeller spectrum 
match the measured values but the high frequency harmonics are 
underpredicted. This underprediction is believed to be due to the omission of 
unsteady loading effects in the predictions. The farfield or flyover noise of 
the subsonic propeller is scattered by atmospheric and ground effects but the 
general trend of the data indicates overprediction of farfield levels. 
Transonic propeller noise measured on the fuselage of the aircraft is 




4 ^ 



2 


significantly influenced by the refraction effects of the boundary layer on 
the fuselage. When these effects are included and when the power is matched, 
the transonic predictions agree with the data except in a small region just 
behind the propeller plane on the aircraft surface. It is believed that 
scattering effects may be the cause of this discrepancy. 


/ 
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INTRODUCTION 

Propeller noise prediction is based on two disciplines: aerodynamics and 

acoustics. NASA has developed a computer system for propeller noise predic- 
tion using the assumption that these disciplines are separable, that is, that 
the flow field can be separated into an aerodynamic part and an acoustic 
part. This separation allows the computations to be made sequentially. 
Classical aerodynamic theory is used to find the surface pressures and 
frictional stresses on the blade surfaces and then acoustic theories are used 
to predict the noise. 

The prediction system is diagrammed in figure 1. The different computa- 
tional tasks are assigned to independent blocks of computer code called 
functional modules. These modules are managed by an executive system, the 
ANOPP executive program (reference 1). One group of modules deals with the 
aerodynamic computations. Not shown are modules which compute atmospheric 
properties and flight dynamics. These modules are described in reference 2. 
The end products of the aerodynamics modules are the propeller blade motions 
and loads. Motions include only the aircraft motion and rotational effects. 
Flexing and vibration, while possibly important, are not included at this 
time. The loads are the pressure and frictional stress on the propeller blade 
surface. These loads are generally a function of both surface position and 
time. Given the blade motions and loads, it is theoretically possible to 
predict the noise. As a practical matter, however, this prediction is impos- 
sible. The loads are really non-stationary random processes and there is 
presently no feasible computational procedure which will produce a complete 
description of the noise. The approximation is made that the noise may be 
divided into two parts: discrete and broadband. The discrete tone noise is 
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computed directly from the blade motions and loads while the broadband noise 
is estimated by semi -empi ri cal methods. The tone noise computations are made 
entirely with time domain techniques while broadband noise estimates use a 
blend of frequency and time domain concepts. 

Noise predictions are needed in both the near and far field. Nearfield 
noise predictions are needed to find the noise transmitted to the interior of 
an aircraft. There are two important nearfield effects on propeller-generated 
noise. These effects are the refraction caused by nonuniforn flow over the 
aircraft and the scattering by the aircraft body. Farfield effects are needed 
to accurately predict community noise. These effects are the atmospheric 
attenuation and refraction, and the ground effects of reflection and absorp- 
tion. Further information on the farfield propagation modules may be found in 
reference 2. Modules for the nearfield effects will be described in more 
detail here, 

Agreement with experiment is the ultimate goal of a prediction system. 

Any prediction method, a guess, a curve fit, or a solution to a partial 
differential equation, is acceptable if it agrees with experimental data 
according to some objective rule. The question is how many experiments and 
how good must the agreement be to prove that a prediction is correct. This 
question raises the second question: prove it to whom? You cannot prove to a 

member of the Flat Earth Society that the earth is round. 

Despite these minor difficulties, we gather data from experiments in the 
hope of proving that the prediction system is correct. These are indicated in 
figure 1 as flight data and tunnel data to denote the two types of facilities 
most frequently used to conduct the experiments. The flight data are usually 
but not always full scale. Results will be shown later for two propellers: 
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one subsonic design and one transonic design. The subsonic propeller is made 
by the Hartzell Company and its noise was measured by a wing-mounted boom 
microphone in flight. The transonic propeller is the NASA SR-3 Prnpfan 
design. This roughly quarter-scale propeller was tested in flight mounted 
atop a Jetstar aircraft. 


SYMBOLS 


radius of perfect circle 
Joukowski transform parameter 
blade section chord 

blade section drag coefficient re section dynamic pressure 
and chord 

energy disipation coefficient 

skin friction coefficient re section dynamic pressure 

blade section lift coefficient re section dynamic pressure 
and chord 

Theodorsen transform coefficients 
coefficient of pressure re section dynamic pressure 
wall shear stress coefficient re local dynamic pressure 
ambient speed of sound 

Prandtl tip vorticity function 

boundary layer shape factors, i ^ j = 1, 2, 3 

Hankel functions of the first and second kinds 

unit imaginary number 
advance ratio Voo/nD 
wave numbers 

loading intensity vector 
empirical broadband correlation length 
tip helical Mach number 
radiation Mach number 
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tip rotational Mach number 
freestream or forward Mach number 

circumferential harmonic number 
blade section Mach number 

normal to blade surface 
acoustic pressure 
cylindrical coordinates 
propeller disk radius 

radiation vector x-y 

unit radiation vector 
surface area 

empirical broadband spectrum 
time 

tangents to blade surface in spanwise and chordwise directions, 
respectively 

period of signal or transfer function 
local Mach number in boundary layer 

velocity vector 

complex velocity field around perfect circle 

distance along chord or distance along airfoil surface measured 
from stagnation point 


observer position in media-fixed reference frame 


y source position in media-fixed reference frame 

2 complex plane of the airfoil section 


Greek Symbols 
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a 

blade section angle of attack 


a' 

angle-of-attack perturbation 


P.75 

blade pitch at 3/4 span 



boundary layer displacement thickness 


&2 

boundary layer momentum thickness 


^3 

boundary layer energy thickness 


e 

polar coordinate angle difference between near circle and 
perfect circle 


fii » 02. na 

pitch axis coordinates fixed to propeller, nz is the pitch 
axis and na is the forward shaft axis 


?• 

complex plane of the airfoil section represented as a near 
ci rcle 


e 

angle between normal vector and radiation vector 


r 

circulation around airfoil section or curve of intersection of 
collapsing sphero and propoller blade surface 


X 

advance ratio M /M^ 
t 


> ^2 

blade surface coordinates. varies with span, varies with 

chord ^ 


P«, 

ambient density 


a 

propeller solidity 


T 

retarded time 


<!> 

blade section inflow angle 



elliptic blade thickness function 


(0 

angular frequency 


n 

propeller angular velocity 
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Subscripts 

h 

i 

L 

r 

s 

S 

tr 

T 


helical 
i nduced 
leading edge 
refraction effect 
scattering effect 
stagnation point 
transition point 
trailing edge 
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AERODYNAMIC PREDICTION 
Geometry 

While geometry is given rather than predicted, it is worthwhile to care- 
fully consider the way geometric data are utilized. Two coordinates are 
required to specify a point on a surface. It is important to choose surface 
coordinates such that the functions of these coordinates will be single-valued 
and free of singularities at least up to their second derivatives, that is, 
functions of the surface coordinates will be of class C^. The selected sur- 
face coordinates should also provide a convenient computational grid. 

Figure 2 shows the basic coordinate system yiy-^y^ used to describe the 
propeller. At time t = 0 , a rotating coordinate system 01020$ is congruent to 
the yiy2y$ system. The shaft axis is o$, the blade pitch axis is 02* snd the 
oi axis completes an orthogonal triad. Sections A-A through the blade at 
constant span positions give the first surface coordinate 



( 1 ) 

( 2 ) 

( 3 ) 

( 4 ) 


is used to introduce the second surface coordinate Kz* The blade surface 
elliptical coordinate t|>(C2) resulting from this transformation is shown in 
figure 3 . The blade surface is described by the function ij){ 52 )* The blade 
surface is unwrapped by the Joukowski transform if the second surface coordi- 
nate ?2 chosen to equal ri2» The surface function t( 5 i» ^z) is single 
valued, continuous and slowly varying in the surface coordinates ^2 3 S 
shown in figure 3 . The computation grid stretches the region near the 
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leading edge, ^2 ~ so that aerodynamic functions such as the coefficient of 
pressure Cp will be slowly varying in this region. In this computational 
space bicubic splines are suitable interpolating functions and will be used in 
all subsequent computations. 


Potential Field 

The potential flow field around each airfoil section = constant is 
given by a conformal transformation of the flow around a perfect cylinder 
(reference 3). The blade geometry analysis has already produced part of the 
desired transformation by mapping the airfoil section in the z-plane as shown 
in figure 3. Recall that the airfoil thickness function t(Ci, ^ 2 ) was gener- 
ated by inverting the Joukowski transformation. 



This transformation may now be used directly to map a given flow around the 
near-circle into a flow around the airfoil section. 

Iheodorsen's transformation (Ref. 4) maps the c-plane of the perfect 
circle into the c '-plane of the near circle 

C' = C exp { I } (6) 

n=l c" 

The constants Cp in Theodorsen's transformation are found from the shape of 
the airfoil in the c '-plane. After numerically solving for these constants, 
it is found that the trailing edge point of the airfoil is displaced by a 
small angle e-;- from the real axis of the c-plane. This point is required to 
be a stagnation point for the flow around the cylinder in order to satisfy the 
Kutta condition that trailing edge velocities are finite. 


The complex flow function around the perfect circle is 


W(c) = M (e-''“ - e^“) + i-L 

^2 2ttc 


( 7 ) 


where M is the local Mach nunber (for the section), a is the angle of attack, 
a is the radius of the circle, and r is the circulation. The trailing edge 
stagnation point is 


le. 


^ST = ^ ® 


which, with W(c^^) = o, gives the circulation r as 
r = 4 ti a M si n(a-e^) 

The leading edge stagnation point, a second solution to the equation W(?^) 

i (2a- e j+tt) 


1 s 


tSL = a e 


( 9 ) 


The cross product of the velocity vector and the circulation vector gives the 
lift 

^ sin(a - e^) ^0) 

The coefficient of pressure is found from the velocity in the z-plane of 
the ai rfoi 1 , 

s ■ > - 1 ^ r (11) 

The complex velocity in the z-plane is found by using the derivatives of the 


transformations 
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At the stagnation points in the ?-p1ane, the limit of equation (12) is used. 

Since the Joukowski transformation has a singularity at the trailing edge, 

W(c) will have a finite non-zero limit at this point. This causes the 
coefficient of pressure at the trailing edge to have magnitude less than one. 

Figure 5 shows the blade section data computed by the blade section 
potential flow analysis. The coefficient of pressure Cp has a maximum of 
unity at the leading edge stagnation point C 2 $» The lower surface Cp is 
generally positive, decreasing to a small value near the trailing edge ?£ ~ 2ir 
where the airfoil surface velocity is near the free-stream velocity. Hie 
upper surface Cp may become negative with high velocities around the highly 
curved leading edge. The upper surface Cp then approaches a small negative 
value at the trailing edge C 2 = 0. The lift coefficient Cl has a slope of 
approximately 2n when plotted as a function of a. 

The function of the blade aerodynamics module is summarized as follows. 
Input is the blade shape function 'J^(Ci,C 2 ) afid the parameters a and M. The 
module maps the flow around a perfect cylinder into the flow around the blade 
using the Theodorsen and Joukowski transforms. Outputs are the section lift 
coefficient and leading edge stagnation point 52 S a function of span 
position Cl and the parameters a and M. The output coefficient of pressure 
depends on surface position (Ci,C2) The parameters a and M. The 
stagnation point location and the coefficient of pressure are used in the 
boundary layer analysis which follows. 


14 


Boundary Layer 

The boundary layer module computes the two-dimensional boundary layer on 
each airfoil section as shown in figure 6. The blade arc length x is measured 
from the stagnation point on the leading edge. The initial portion of the 
boundary layer is laminar. Transition to turbulent flow occurs near the point 
where the external velocity U(x) is a maximum. This turbulent layer continues 
to the trailing edge unless separation occurs. 

The governing equations for the boundary layer thicknesses are the 
integral momentum equation for 62(x) and the integral energy equation for 
63 (x). The displacement thickness 61 (x) is related to the momentum and energy 
thicknesses through the assumed profile U(y) for the boundary layer. 

The governing equations for 62 (x) and 63 (x) are 


(x ) 

dx 


+ ( 


u 


dU. 


62 (x) = 


(14) 


+ (-g) 63 (x) = Cp (15) 

The shape factor H12 is a given function of the thicknesses 62 and 63 and the 
external velocity gradient The coefficient is the local wall shear 

stress coefficient and the coefficient Cg is the energy dissipation 
coefficient. 

In the laminar layer, Cj is known from the assumed boundary layer 
velocity profile and equation (14) can be integrated to find 62(x) without 
solving for the energy thickness. Holstein and Bohlen's method (reference 5, 
Chapter X) is used to integrate this equation. Transition is assumed to occur 
when the external velocity is a maximum, that is, where 
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dlUx^-p) 

~dx* " ” (Ifi) 

Both equations must be integrated in the turbulent region. Truckenbrodt ' s 
method (reference 5, Chapter XXII) is used for this purpose. The coefficients 
Cy and Cp are given by empirical functions in the turbulent region. 

The section profile drag coefficient is computed by Young and Squire's 
method (reference 5, Chapter XXV). In this method, the wake thickness at 
infinity is estimated by the empirical formula 

IJ(Xt) 

*^2^ = >S2(Xy) ( ji j (17) 

00 

and the section drag coefficient, referred to the chord, is 

^d " (18) 

The boundary layer module computes the skin friction coefficient (see 
figure 7) Cy for use later in the blade loading module. It provides the 
drag coefficient C(j for propeller performance analysis and for the computa- 
tion of lifting line (compact source) loads. The trailing edge thicknesses 

(Xy) 3 nd 62 (Xy) are used in scaling laws for trailing edge broadband 
noise. 


Propeller Performance 

Computation of the propeller performance depends on a solution for the 
induced velocity field. This induced field at a blade section is shown in 
figure 8 . Propeller performance is predicted by Lock's method (reference 6 ). 
The blade section aerodynamic module gives tables of section lift functions 
Ct(Ci, 01 , M). The boundary layer module gives tables of Cd(Ci, a, M). 
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These tables are used in the prediction of propeller performance coefficients 
Cp and Cj for given tip and forward Mach numbers Mt and M of the 
propeller. 

The helical Mach number at a blade section (see figure 9) Mf, is known 
from its components. These components are the axial Mach number of the pro- 
peller and the rotational Mach number at the section ^jMt, where Mt 
is the tip rotational Mach number Rn/Coo. The local induced Mach vector M-f 
must be found to compute the total local Mach vector by 

^ ^ ^ ( 1 ?) 
This vector is represented by its magnitude M and the inflow direction angle 
^ as shown in figure 9. Differential components of lift dL and drag dD on the 
section are rotated through the inflow angle ip to give the differential thrust 
force dT and torque force dQ. These differentials are then integrated over 
the blade length and converted to performance coefficients; the torque or 
power coefficient Cp and the thrust coefficient Cj. The two components of 
the induced Mach vector are supplied by solving the two equations for the 
change of momentum in the far wake of the propeller. The increase in axial 
momentun, the value downstream minus the value upstream, is equal to the 
thrust force. The increase in angular momentum is equal to the propeller 
torque. These balance equations are 

2Msin4. (Msin4>-M^) F(Ci,X) = | M^aC^^ )[C^cos<),- C^sinfl 

2Msin4. {CiM^-Mcos.},) F(5i,X) = M2o(5j[C^sin((.+CjCos<|.] 


( 21 ) 
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The factor of 2 on the left of these equations represents the fact that 
velocities in the far wake are twice those at propeller disk. The factor M 
sin(|) represents the mass flow through the disk. The factor (Msin<{) - M^) in 
the thrust equation is the axial induced Mach number and the factor (?iMt - 
McoS(|)) is the angular induced Mach number. The Prandtl circulation function 
(reference 7) F(Ci,X) is derived from vorticity theory. It is defined by 

F = I Arccos { exp | (1 - ] (22) 

and is an approximate way of representing Goldstein's circulation function 
(reference 8). The above solution procedure was developed in 1930 by C. N. 

H. Lock (reference 6). 


Loading 

Blade loads are computed by combining results of the aerodynamics, 
boundary layer, and performance modules as shown in figure 10. The lift and 
drag coefficients are three-dimensional tables in terms of span position 5^, 
angle of attack a, and Mach number M. The surface stress coefficients Cp 
and Cf are four dimensional tables in terms of surface position and ^ 2 * 
angle of attack a, and section Mach number M. The performance analysis 
produces actual values a(Ci) and M(?i) for the propeller operating conditions 
"Oiere may be an additional angle-of-attack perturbation a'(5i,t) 
due to small nonuniformities in the propeller inflow. 

When the stress and loading coefficient tables are interpolated with 
these functions, the coefficients became actual time-dependent values 

a(?i,t), M(Ci)] (23) 


Cd(^i.t) = C^C^i, a(?i,t), M(Ci)] 


(24) 


Cp(5i.C2,t) = CpC^i, Kz, cx(5i,t), M(5i)] 
Cf(5i,52,t) = C^CCi, ?2. a{5l.t), M(5i)] 

where 
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( 26 ) 


Replacing the parametric arguments a and M by actual arguments reduces the 
lift and drag tables to two dimensions and the stress coefficient tables to 
three dimensions. If steady loadings are assumed, the loading tables are 
reduced by a further dimension. These loading tables are passed to the 
discrete noise prediction module for noise prediction. The line loads 
3nd C(j{Ci,t) are used in compact source theories and the 
distributed loads Cp(5i,^2>^) 3nd Cf(C 2 >C 2 >t) are used in the general 
non-compact source theories. A similar transformation of the boundary layer 
thicknesses is used to supply data to the broadband noise prediction module. 
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ACOUSTIC PREDICTION 
Di screte Tone Noise 

Subsonic Noise .- Discrete tone noise is predicted by Farassat's method 
(reference 9). Farassat's equation for the noise of a subsonic propeller is 

M +£ X, 


4Trp(x,t) = 


. at 


/ [ ds + / [ 1 dS 

blade t blade 1-M^ 


(28) 


The terms in equation (28) are illustrated in figure 11. Pressure is computed 
at a particular observer position x and time t. Two integrals are evaluated 
to give the total pressure. The integrals are over the surface area of the 
propeller blade and the integrands are evaluated at the time x when the sound 


is emitted at the surface position y(x). The radiation vector r is the 
difference between the observer and source positions 


r = x(t) - y(x) (29) 

The normal vector n and surface area dS are given by 

ndS = 1^ X If- (30) 

0^1 0^2 

and the velocity vector is 

V = s Xn (31) 

The loading vector is 

1 = ^ (± c /2 + Cpn) (32) 

where the positive sign is used on the lower surface between the stagnation 
point on the lower surface and the trailing edge. The base vector ^2 is 

tangent to the surface and in the chordwise direction. A unit vector in the 

t 

radiation direction r is used to define 



t 

r • V 
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' = * • * (34) 

Integration of equation (28) is straightforward except for the solution 

for retarded time. The basic retarded time equation requires the distance 
from the source to the observer to be compatible with the propagation time 

r . ? = j^(t) - ^( x )|2 = (t _ x )2 ( 35 J 

In the case of a propeller moving along its own axis, the retarded time equa- 
tion can be reduced to 


+ 2BQ(t-x) + Cq + Cl cosx = 0 (35 

Equation (36) has the appearance of a quadratic equation in (t-x) except for 
the coefficient C = Cq + Cicosx. It can be shown, however, that there is a 
single real solution x < t to this equation as long as the motion of the 
propeller is subsonic. 

Transonic noise.- Supersonic noise is computed by Farassat's collapsing 
sphere method. The subsonic equation cannot be used on any portion of the 
blade when may exceed unity because of the |l-Mr| singularity. In 
addition, the retarded time equation has multiple roots. The collapsing 
sphere method is illustrated in figure 12. The collapsing sphere intersects 
the blade surface in a curve called the r-curve. Farassat has shown that 



( 37 ) 
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The integrals for noise are thus evaluated by choosing a set of times x and 
integrating first along the r curves with fixed source time and then over 
source time. This method was developed by ^^ystrom and Farassat in reference 
10 . 

Recent advances . -The time derivative in Farassat's acoustic equation 
increases the computation time because at least two integrals must be 
evaluated to numerically compute the derivative. Numerical differentiation 
also introduces some spurious wiggles in the pressure signature which appear 
as increases in the higher harmonics of the transformed signal. Recently, 
Farassat (reference 11) has taken the derivatives inside the integral for both 
the subsonic and supersonic cases. The differentiation under the integral has 
shown that the noise depends on blade surface curvatures. The full 
implications of this exciting new result are not completely known at this time 
and are a subject of continuing research. 
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Broadband Noise 


Broadband noise is generated by turbulence convected past the trailing 
edge of the airfoil as shown in figure 13. Hie mean-square pressure spectrum 
is given by Schlinker and Aniet as 


<p2(J.„)>= / 1 |l|2, 5(„,J 

trail inn 9-n ' ' j 


where 


trailing 

edge 


2tt 


X U 


a = [x^ + B^(y^ + z2)]V2 

The function L is the "effective lift" function derived by Aniet (reference 
12) for an observer at retarded position (x, y, z). The function £y(o)) is 
an empirical correlation length 




2.1 U 


where is the turbulence convection velocity which is about 0.8 U. 

The spectrum function S(oj) is an empirical function for the blade surface 
pressure spectrum at the trailing edge. 

S(“) = 2xl0"5 (i F((o) (42) 


F(o)) = 33.280) (1 - 5.49o) + 36. 7o)^ + O.lBlo)*^)"! 


23 


ftearfield Effects 

Two effects are present in the nearfield of the propeller which may 
significantly alter the received noise. These effects are the scattering by 
the aircraft’s wings and fuselage and the refraction by the boundary layer on 
the surface of the aircraft. 

Scattering. The scattering effect is illustrated in figure 14 for a 
cylindrical fuselage. The free field levels on the surface of a cylindrical 
fuselage are calculated by one of the previously described methods for 
discrete noise. This incident field pi(o),0,x) is transformed by 

2tt “ -i (k X + me) 

P^(w,m,k^) = -|^ / / e p. ((o.e.x) dxde (44) 

0 

to give the incident field in a wavenumber space. The solution for the total 
pressure on the fuselage surface can then be found by superimposing the 
general solutions for incident and scattered cylindrical waves such that the 
boundry condition 

3p(a),m.k ) 

— 3- = 0 (45) 

'r=a 

is satisfied. The result for the total surface pressure can be given by a 
transfer function. 

p^(o),m,k^) = 2T^(a),m,k^) p. (u.m.k^) 

where 

iTTk a (1)’ (2) -1 

T 3 (».m,k^) = 1--^ (k^a) H„, (k,a)} (47) 


and 




( 48 ) 
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Tlie factor of 2 in equation (46) is shown explicitly to represent the 
effect of pressure doubling. With this form, the transfer function for scat- 
tering approaches unity for high frequencies as may be seen from the 
asymptotic forms for the Hankel functions. Following the computation of total 
pressure pt in wavenumber space, an inverse transform is used to find the 
surface pressure as a function of position (0,x). 

Refraction .- The boundary layer velocity profile alters the sound pres- 
sure level on the surface of the fuselage by turning the waves propagated 
upstream away from the surface and turning waves propagated downstream into 
the surface. This effect is depicted in figure 15. This refraction effect is 
small for aircraft in low speed flight but becomes significant at the higher 
subsonic Mach numbers. 


A simple model of the refraction effect neglects scattering and uses the 
two-dimensional wave equation in a sheared flow as the basis for finding a 
transfer function. The sheared-flow wave equation is 

(a,_Uk^) + 2(^) l<x ^ + (“-Uk^)[(o)-Uk )2_ k^] p = 0 (49) 

dy^ 


This equation can be integrated, given the boundary layer velocity profile 
U(y), from the surface where y = 0 to the edge of the boundary layer where 
y = 6. Initial conditions at the surface are an assumed unit pressure 


and 


P(0) = 1 



(50a) 

(50b) 


Integrating equation (49) with initial conditions (50) gives the pressure p^ 
and velocity vg at the edge of the boundary, layer. These results must be 
scaled to match the pressure in the known incident wave field. The process of 
matching the incident field produces a transfer function for refraction 
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T^(o),k,k^) = Pg(o),kJ - (oj-lik^) Vg((o,k^) (51) 

Numerical integration of the sheared flow equation breaks down at the regular 
singular point (u - Uk^) = 0 which corresponds to the coincidence of wave 
speed with flow speed in the boundary layer. Special techniques beyond the 
scope of this paper have been used to integrate the sheared-flow wave equation 
(49) in these cases. 

Combined effects .- The combined effects of scattering and refraction may 
be found by integrating the sheared-flow equation in cylindrical coordinates 
to find the cylindrical wave pressure pg and velocity v^ at the edge of 
the boundary layer. Matching the external field then gives a combined 
transfer function 


T „(o),m,k ) = 
sr' ’ * x' 




(52) 


z = k^(a+6) 


where 


(53) 
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Farfield Effects 

Noise propagated to the farfield is modified by the effects of atmos- 
pheric attenuation and ground reflections and attenuation. Since phase infor- 
mation tends to be lost after propagation over long distances, farfield noise 
effects are applied to the mean-squared pressure spectrum of the noise rather 
than to the pressure itself (reference 2). 

Atmospheric attenuation .- The effect of atmospheric attenuation is 
computed by the ANSI method described in reference 2. The mean-squared pres- 
sure on a source sphere of fixed radius is reduced by the spherical 
spreading effect and atmospheric attenuation. The transfer function for 
attenuation Tg is a decaying exponential function exp(-2ar). The attenua- 
tion rate a is a function of frequency. The lower frequencies are dominated 
by the effect of nitrogen relaxation. Mid-frequencies are usually dominated 
by 0 )^gen relaxation effects and the higher frequencies are dominated by the 
classical absorption effects of conductivity, viscosity, and rotational modes 
of molecular vibration. The ANSI-proposed standard method is used for the 
calculation of attenuation effects for both standard and non-standard 
conditions. 

Ground effects .- The farfield noise is reflected and attenuated by the 
ground. This effect may be represented in the transfer function Tg developed 
by Pao, Wenzel, and Oncley in reference 13. This transfer function is based 
on the complex reflection coefficient Re'i't’ for a spherical wave over an 
impedance plane. The factor R is the magnitude of the reflection which is 
equal to or less than unity. The factor ({> is the phase shift between the 
reflected and incident waves. The magnitude of the noise at the observer 
depends on the difference in the lengths of the direct ray path and the 
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reflected or Image ray path. The magnitude of the received noise is 
diminished slightly by the loss of coherence of the direct and reflected 
signals. This coherence loss also depends on the path length difference. 
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COMPARISONS TO EXPERIMENTS 
Subsonic Propeller 

Noise predicted by this system has been compared to noise measured during 
the flight of the twin-engined aircraft shown in figure 16. The aircraft was 
fitted with a wing-mounted microphone boom which could be moved to measure the 
noise in front of and behind the propeller plane. Noise was also measured on 
the ground with level flights over the measurement point. 

Propeller performance .- An intermediate check of the prediction system 
was made by Block and Martin (reference 14) using a 1/4-scale model of the 
twin aircraft's propeller on the propeller test stand (figure 17) at Langley 
Research Center. 

Computed and measured power coefficients are shown in figure 18 as a 
function of advance ratio J = Voo/nD and 3/4 span pitch setting B. 75 . The 
predicted power coefficient is near the measured values at the lower pitch 
settings but rises above the measured power at the higher settings where it is 
believed that the propeller is in a condition of partial stall. Similarly, 
the predicted thrust coefficient Cy (figure 19) is above the measured value 
with fair agreement at the low pitch settings and poor agreement at the high 
pitch and thrust values. 

Effective pitch .- Pitch was not measured during the flight of the air- 
craft when the noise was measured. In order to find an effective value for 
the pitch setting B.75 to use in noise predictions, the measured power in 
flight was used to compute the power coefficient Cp. An effective pitch 
setting (B,7s)eff was then found such that the predicted power coefficient 
matched the measured value. This effective pitch will be less than the pitch 
setting on the 1/4-scale model for the same Cp. 
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Propeller noise . Nearfield noise measured in flight is shown in figure 
20. The measured noise data were sampled using data blocks of 96 time points 
per shaft revolution which gave 32 points per blade for the three-bladed 
prop. Fifty blocks of data were ensemble averaged to find the mean signal 
p(t) shown in figure 20. Signal number 1 is for one-third of a propeller 
revolution. Signals 2 and 3 which complete the revolution are similar. The 
standard deviation o(t) was nearly constant, indicating that the measured 
noise could be decomposed into discrete and random parts 


p(t) = p(t) + p'(t) (54) 

where p' (t) is a stationary random signal. Sound pressure spectra were gene- 
rated for each data block and for the ensemble-averaged data. The 

discrete spectrum jp(a))j^ and the random spectrum jp'(w)j^ are added to 
give the total spectrum jp{tD)j^. 

|p{o))|2 = |p{o))|2 + |p'(o))|2 (55) 

The spectral data shown in figure 21 are relative to the overall mean-squared 
pressure 


<p^> = 1 im 

T^oo 



p^(t)dt 


(56) 


The overall mean-squared pressure in figure 21 is within 1.8 dB of the 
measured value. This agreement may be due in part to the matching of power 
through the effective pitch setting procedure. 

The predicted spectrum agrees well with the data for the first few 
harmonics. Starting at about the 5th harmonic, the measured discrete spectrum 
shows a cyclic pattern suggestive of cancellation and reinforcement which is 
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not predicted. This pattern gives levels from the 10th and 15th harmonics 
which are as much as 20 dB above the predicted values. A significant part of 
this underprediction is probably due to the omission of unsteady loads in the 
analysis. Despite these high-frequency errors, the predicted A-weighted sound 
pressure level will be in good agreement with the measured values because of 
the rapid decay of the levels with frequency. 

The error in the discrete spectrum may be due to any of several effects, 
including the obvious possibility of an error in the prediction method. Scat- 
tering from the wing and fuselage was not included in these computations. 
Computations are made for one propeller and there may be a small contribution 
from the second propeller which contributes to the error. 

In making the prerdictions, it was assumed that all blades are identical, 
whereas, in practice, there are slight blade-to-blade differences which may 
result in the reinforcement patterns seen here. There is regretably no way to 
examine these possible sources of error without further detailed 
experimentation. 

Farfield noise from the flyover of the twin prop aircraft is shown in 
figures 22 and 23. The limited amount of data shown in figure 22 are for 
flyovers with the reception angle 6 = 90“ being the directivity angle when the 
sound is received. The prediction is 5 to 6 dBA above the average of the 
measured levels. The reason for the overprediction is shown in figure 23. 

The leveT for the first harmonic, which dominates the predicted source spec- 
trum, is indicated by the data to be less than the second harmonic. The 
reason for this effect is unknown and a much larger data set should be 
examined before any conclusions are drawn. Flyover noise data for propellers 
typically scatter over a range of about 10 dB (see reference 15) so that 
statistical methods must be used to assess the accuracy of predictions. 
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Advanced Turboprop 

Predictions have been made for the eight-bladed advanced turboprop 
propeller designated as SR-3. Since this propeller was operating at transonic 
conditions, the elementary aerodynamics and loading methods described earlier 
could not be used. A locding dictribut icn by Ihe propeller 

manufacturer and this distribution was adjusted such that the computed power 
matched the power measured during the flight test of the propeller. 

The propeller was tested in flight atop the Jetstar aircraft shown in 
figure 24. Performance data were also obtained in the wind tunnel as shown in 


figure 25. Microphones were mounted in the surface of the fuselage along a 
line under the axis of the propeller. Measured noise data were reduced using 
ensemble-averaging techniques as described earlier. Only the results of the 
discrete noise data and predictions will be shown here. 


Figure 26 shows a comparison of the measured and predicted noise at the 
cruise Mach number of 0.80. The helical Mach number of the blade tip is 1.13 
for this example. Levels shown are for the blade passing harmonic which, 
remembering the subsonic results, is expected to give the best agreement with 
the measured data, tiystrom and Farassat's PROPFAN (reference 10) program was 
used to predict the freefield noise at the top of the boundary layer. The 
boundary layer profile was measured with rakes mounted on the fuselage of the 
Jetstar. The two-dimensional transfer function for the boundary layer was 
then used to predict the noise on the surface beneath the boundary layer. It 
is apparent that the boundary layer has a significant effect on the surface 
noise at this flight speed. The predictions agree better with the data both 
in front of and behind the plane of the propeller. 

The causes of the discrepancy between predicted and measured levels on 
the Jetstar just behind the prop plane are being investigated at this time. 
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The scattering effect given by equation (47) will be evaluated first to see 
how scattering modifies the predicted surface pressures in the absence of a 
boundary layer. Following this the combined effects of scattering and 
refraction will be computed using the transfer function (52). 

CONCLUDING REMARKS 

NASA's ANOPP system for predicting propeller noise has been described 
here. The objective of this system is to predict noise directly from the 
shape and motion of the propeller blades. This objective requires aerodynamic 
computations to be made as a basis for the acoustic computations. The aerody- 
namic theories used here are simple classical methods. 

These methods overpredict the performance of the propellers which results 
in an error in the acoustic predictions. An empirical correction procedure of 
matching predicted and measured power was used here to compensate for the 
known error in aerodynamic performance. The acoustic predictions give good 
agreement with the low-frequency discrete noise when they are adjusted by the 
power matching procedure. It is believed that improved aerodynamic theories 
will give good low-frequency predictions for subsonic propellers without an 
empirical correction for power. 

The accuracy of the acoustic predictions decreases with increasing 
frequency. Errors of 10 dB or more may occur at or above the lOth harmonic of 
the blade passing frequency in the case of subsonic propellers. When the 
propeller has nearly steady loading, these errors do not seriously affect 
integrated measures of noise such as the A-weighted sound pressure level 
because the A-level is dominated by lower harmonics. The prediction methods 
used here are believed to be adequate for subsonic propellers operating in a 
tractor configuration where the loads are nearly steady. There is a need for 


33 


further investigation of noise from subsonic propellers operating as pushers 
where the unsteady loads will increase the levels of the higher harmonics 
relative to the fundamental. 

Nearfield effects of scattering and boundary layer refraction 
significantly alter the noise levels of the fuselage of the aircraft, llie 
boundary layer refraction effect reduces the surface noise in front of the 
propeller and increases the noise behind the propeller. This effect is 
significant for an aircraft with high subsonic Mach number such as the n.80 
Mach nunber envisioned for advanced turboprop-powered aircraft. The 
two-dimensional model of the boundary layer fails to explain a significant 
discrepancy between predicted and measured data just behind the plane of the 
SR-3 turboprop on the Jetstar aircraft. Three-dimensional refraction and 
scattering effects are being studied in an attempt to resolve this 
discrepancy. 

Farfield effects of atmospheric attenuation and ground reflections modify 
the noise measured during the flyover of a propeller aircraft. The largest 
effect is that of ground reflection which introduces the "ground-dip" in the 
measured flyover noise. Available theories and methods are adequate for 
predicting these effects when the ground reflection angle is large. Although 
measurements from outdoor flyover tests have a typically large scatter of 
data, the flyover noise from propellers can usually be predicted within a 
standard deviation of about 4 dB. 
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Figure 1.- ANOPP propeller noise prediction system. 
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Figure 4.- Potential flow mapping for blade section. 
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Figure 5.- Propeller blade lift and pressure coefficient functions. 
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Figure 6.- Airfoil boundary layer analysis. 
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Figure 7.- Propeller blade drag and skin friction coefficient functions. 
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Figure 10.- Blade lift and drag loading function. 
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Figure 11.- Farassat's method for predicting discrete subsonic propeller noise. 
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Figure 13.- Schlinker and Aniet's method for predicting trailing edge broadband noise. 





Figure 14.- Scattering of nearfield source noise by a cylindrical fuselage. 
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Figure 15.- Fuselage boundary layer refraction of nearfield source noise. 




Figure 17. Quarter-scale subsonic propeller in test chamber. 
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Figure 20.- Inflight nearfield subsonic propeller noise signal. 
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Figure 21.- Inflight nearfield subsonic propeller noise spectrum. 
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Figure 22.- Twin subsonic propeller aircraft flyover noise levels. 
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Figure 23.- Twin subsonic propeller aircraft flyover noise spectrum. 








Figure 24.- Jetstar/Propfan flight tests 




1 « ♦ 



I 


f 



1. Report No. 2. Government Accession No. 

NASA TM 85636 

3. Recipient's Catalog No. 

4. Title and Subtitle 

Propeller Noise Prediction 

5. Report Date 
May 1983 

6. Performing Organization Code 

505-40-12-01 

7. Author(s) 

william E. Zortomski 

8. Performing Organization Report No. 


10. Work Unit No. 

9. Performing Organization Name and Address 

NASA Langley Research Center 
Hampton, VA 23665 

11. Contract or Grant No. 

13. Type of Report and Period Covered 

Technical Memorandum 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, DC 20546 

14. Sponsoring Agency Code 

15. Supplementary Notes 


Presented at Acoustical Society of America Meeting, May 10, 1983 in 
Cincinnati, Ohio. 


16. Abstract 

Analytic propeller noise prediction involves a sequence of computations culminating 
in the application of acoustic equations. This paper describes the prediction sequence 
currently used by NASA in its ANOPP (Aircraft Noise Prediction) program. No attempt 
is made here to review the state of the art of noise prediction. Some elements of 
this sequence represent classic results while others represent the most recent publi- 
cations. The elements of the sequence are called program modules. The first group 
of modules analyzes the propeller geometry, the aerodynamics, including both potential 
and boundary layer flow, the propeller performance, and the surface loading distribu- 
tion. This group of modules is based entirely on aerodynamic strip theory. The next 
group of modules deals with the actual noise prediction, based on data from the 
first group. Deterministic predictions of periodic thickness and loading noise are 
made using Farassat's time-domain methods. Broadband noise is predicted by the 
semi-empirical Schl inker-Amiet method. Near-field predictions of fuselage surface 
pressures include the effects of boundary layer refraction and (for a cylinder) 
scattering. Far-field predictions include atmospheric and ground effects. Comparisons 
are made to experimental data from subsonic and transonic propellers and NASA's future 
directions in propeller noise technology development are indicated. 


17. Key Words (Suggested by Author(s)) 

Propeller perfonnance 
Propeller noise 
Noise prediction 


18. Distribution Statement 

Unclassified - Unlimited 


Subject Category 71 


19. Security Oassif. (of this report) 
Unclassi f ied 


20. Security Classif. (of this page) 
Unclass i f ied 


21. No. of Pages 
61 



*Fof sale by the National Technical Inforination Service, Springfield, Virginia 22161 





















